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Abstract. We consider nonholonomic systems with collisions and pro- 
pose a concept of weak solutions to Lagrange-d'Alembert equations. In 
the light of this concept we describe dynamics of the collisions. Sev- 
eral applications have been investigated. Particularly the collision of 
rotating ball and the rough floor has been considered. 



1. The Description of the Problem 

Let us start from the following model example. There is a solid ball B of 
radius r and of mass m and let its centre of mass coincide with the geometric 
centre S. The moment of inertia relative to any axis passing trough the point 
S is equal to J. 

Give an informal description of the problem. Being undergone with some 
potential forces the ball rolls on the floor and sometimes it collides with a 
vertical wall. After the collision it jumps aside the wall. The wall and the 
floor are rough: the ball can not slide on the floor and along the wall. 

We wish to construct a theory of such a motion in the Lagrangian frame. 
Particularly, we wish to give sense to the term " elastic collision" in nonholo- 
nomic context. 

In physical space introduce a Cartesian coordinate system Oxyz. Let 
(%S, ys, zs) be the coordinates of the point S. 

Suppose that the plane Oxy is a solid and rough floor and the plane Ozy 
is a solid and rough wall. For all the time t > we have z$ = r, x$ > r. 

By C £ B denote the contact point of the ball and the floor. The ball 
can not slide on the floor: 

v c = vs + P,SC] = 0, (1.1) 
here vc is the velocity of the point C, and W is the angle velocity of the ball. 
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When the ball reaches the wall (say by its point G € B) then we also 
have 

v G = v s + P,SG] = 0, (1.2) 
The configuration manifold of the system is M = R 2 x SO(3), where 
(xs,ys) £ R 2 an d an element of 50(3) determines the orientation of the 
ball. We use the Euler angles for the local coordinates in SO(3) . 
Consequently the position of the ball is determined by the vector 

x = (xs,ys, -f,9,^) T - 

Why do we write cp with negative sign will be clear below. 

The wall is a 4-dimensional manifold N = {x$ = r} C M. 

Thus the general construction is as follows. We have a smooth con- 
figuration manifold M, dimM = m and a smooth submanifold N C 
M, dim N = m — 1 (the wall). Both manifolds carry the distributions. 

In the example under consideration the manifold M carries the nonholo- 
nomic constraint given by (jl.ip and the manifold N carries the constraint 
given by (flTTj) and (fL~2|) . 

Let 

x = (xi, . . . ,x m ) T £ M 

be local coordinates in M. 

To determine the distribution at each point x € M introduce a linear 
operator 

A(x) : T X M — > R m ~ l , dmiimA(x) =m-l, xeM 

and the mapping x \— > A{x) is smooth. The subspaces ker^4(x) C T X M 
define an /—dimensional distribution in M. 

To define an s— dimensional distribution in N introduce a linear operator 

B(x) : T X M -> R m ~ s , dimim A(x) = m - s, x £ N. 

The distribution on N consists of the subspaces kerS(x) C T X N. The 
operator B{x) is also a smooth function of x. 

The operators A, B are not uniquely defined: the same distributions can 
be generated by the different operators A, B but we use them because they 
naturally arise in the applications. 

Assume also that 

ker B(x) C ker^4(a;) 

for each x € N. 

The dynamics of the system is described by the smooth Lagrangian L(x, x). 

From the configuration manifold's geometry viewpoint the collisions of 
the rigid bodies was considered in [3]. 

One of the results of this article is as follows. The manifold M is endowed 
with the Riemann metric generated by the kinetic energy of the system. The 
evolution of the system is expressed by the function x(t) £ M. When the 
point with coordinates x(t) collides with the submanifold iV i.e. x(t) € iV it 
then jumps aside obeying the law "the angle of incidence is equal to the angle 
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of reflection" . The authors obtained this law from a system with non-solid 
constraint (wall) by means of some natural limit process. 

These results have been obtained in the absence of nonholonomic con- 
straints. We generalize them to the nonholonomic case. 

In Section 15.21 we consider a rough ball colliding with a floor and obtain 
fromulas which particularly describe the following effect [2]: "A perfectly 
rough ball which conserves kinetic energy behaves in such an unexpected 
way that it is difficult to pick up after it has bounced twice upon the floor, 
and, more bizarre, it returns to the hand on being thrown to the floor in 
such a way that it bounces from the underside of a table." 

2. The Weak Solutions to the Lagrange-d'Alembert Equation 
In the absence of unilateral constraint N a smooth function 
x(t) = (xi, . . . , x m ) T (t) G M, Xi (t) G C 2 [t u t 2 ] 
is the motion of the system if and only if for any function 

supp^fc C (h,t 2 ), ip(t) G kevA(x(t)) (2.1) 
it satisfies the Lagrange-d'Alembert equation 

(^( x (t),x(t))-j t ^(x(t),x(t)))m = 0, te[h,t 2 ] (2.2) 
and the equation of constraint 

x(t) G ker A(x(t)). (2.3) 

Definition 1. We shall say that a function x(t) G H l \t\,t 2 ] is a weak 
solution to the system of Lagrange-d'Alembert equations and the equations 
of constraint iff the equation 

/* 2 /f)T r)T \ 

(— (x{t),x{t))m + QT^mm) * = o (2.4) 

holds for any tp that satisfies h2.1\) and equation 12. 3\) holds for almost all 
t G [h,t 2 \. 

Note that by the Sobolev embedding theorem, the space H l [tQ,t{\ belongs 
toC[t ,ti}. 

In case of smooth function x(t), equations (|2.4|) . (12. 3p are equivalent 
to equations (12. 2j) . (|2.3p . This follows from integration by parts and the 
Lagrange-d'Alembert principle [!],[!]. 

If the motion x(t) contains collisions it is piece- wise differentiable: at the 
moment of collision its first derivative is not continuous and the second one 
does not exist. 

Equations (12. 4h do not contain the second derivative of x(t). Therefore 
the concept of weak solutions is a proper tool to describe the motion with 
collisions. 
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Let us turn to the details. 

Consider a solution x(t) that collides the wall at the moment r G 
i.e. x(t) G N. Correspondingly, one must put 

ip(r) G kei B(x(t)). (2.5) 

We suppose that x(t) G C[ii,i2] and 



x(t) 



x~{t), t€[tl,T], 

x+(t), t€(r,t 2 ] 



and x~(t) G C 2 [ti , r] , x + (i) G C 2 (r, t 2 ] . 
By definition put 

x + (r) = lim x + (r) = lim x + (t). 

We assume that these limits exist. 

The solution x(t) must obey nonholonomic constraint that is 

x ± (t) G ker A(x ± (t)). 

If x~(t) = x + (t) then the derivative x(t) is defined and 

x(j) G ker B(x(t)). (2.6) 

2.1. The Equations of Collision. Introduce the notation 

v ± = x ± (t) G ker A(x{t)). 

Splitting the integral (|2.4p in sum + f^ 2 and integrating them by parts 
obtain 

(^( x (t),x(t))-^(x(t), X (t)))mdt = 0, (2.7) 

(^( x (*)'*(*)) " ^ = 0, (2-8) 

(^(x(r), V + ) - — (x(r),«-))^(r) = 0. (2.9) 

Equations (|2.7p . f|2.8j) express that the functions sc (i) satisfy to the Lagrange- 
d'Alembert equations. By the assumption they also satisfy the equations of 
constraint: x (t) G ker A(x (t)). That is before and after the collision the 
system obeys to the Lagrange-d'Alembert equations and the equations of 
constraint. 

Particularly, if the system is holonomic outside iV (i.e. A(x) = 0) then in 
their domains the functions x^(t) satisfy the Lagrange equations 

^(x ± (t),i ± (t))-|g(x ± (t),i ± (t)) = 0. 



Equation (|2.9p describes the behaviour of the system at the moment of 
collision. This equation is of main importance for us. 
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3. Lemma from Vector Algebra 

The following lemma is mainly used in Section [5j But we placed it here 
because it provides an introduction to the geometry of the next Section. 

Lemma 1. Let L = R m be a Euclidean vector space with scalar product 
given by its Gramian matrix G. And let B be the matrix of a linear operator 
(we denote operators and their matrices by the same letters) 

B:L^R rn ~ s , rtrngB = m- s. 

Let 

L = ker B © W, W _L ker B 

be the orthogonal decomposition of the space. 

Then the square matrix of orthogonal projector P : L — > L, P(L) = W 

is 

P = G- 1 B T (BG~ 1 B T y 1 B. (3.1) 

If an operator 

A : L -> R k 
is such that ker B C ker A then one has 

AP = A. (3.2) 
Particularly, this implies that P(ker A) C ker A. 

Proof. To obtain formula (|3. 1 1) fix an arbitrary vector i£L and introduce 
a linear function /(£) = [Px) T G^. It is clear 

ker B C ker /. 

This implies that there is an operator A : R m_s — > R such that (Px) T G = XB 
and Px = G~ l B T \ T . It remains to find X T from the equation B{x — Px) = 0. 
To obtain formula (|3.2p note that there exists an operator 

such that A = jB. Consequently, formula (j3.2|) follows from (|3. 1 1) . 
The Lemma is proved. 

4. The Natural Lagrangian System 
To proceed with our analysis put 

L = T(x,x,x) — V(x). 

The form 

T(x,£,rj) = ^ T G{x)r], £ = . . . , £ m ) T , rj = (r/i, . . . , r? m ) T 

is the kinetic energy of the system, the matrix G{x) = G T (x) is positive 
definite. It defines a Riemann metric in M. The potential energy V is a 
smooth function in M. 
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By (|2.5j) equation (|2.9f) is reduced to 

T(x(t),v + - v~,u) = 0, u G ker B(x(t)). (4.1) 

From formula (|4.1|) it follows that the difference v + — v~ is perpendicular to 
ker B(x(t)). 

Hypothesis 1. The vector v + depends on the vector v~ by means of a linear 
operator 

v + = R(x(t))v~, R(x(t)) : ker A(x(t)) -> ker A(x(t)). 

Hypothesis 2. The energy is conserved during the collision: 

T(x(t), v + , v + ) = T(x(r), v~ , v~). 

Hypothesis 3. The system is reversible: if x(t) is a motion of the system 
then x(—t) is also a motion. For the collision this implies that: 

v~ = R(x(t))v + . 

The third Hypothesis implies (i?(x(r))) 2 = I. 

Note that if dimker B{x{t)) = dim ker A(x(t)) — 1 then the last Hypoth- 
esis is fulfilled automatically. 

It is reasonable to consider the following decomposition 

T x(r) M = keiB(x(T)) W(x(t)), 

here W(x(t)) is the orthogonal complement for ker B(x(t)), and let 

P : T X(T) M -> W(x(t)) 

is the orthogonal projection. 

Introduce notations Pv = v±, (I — P)v = G ker B(x(t)) and the 
norm |^| 2 = T(x(t),£,£). Then write 

v ± = v^ + 

Theorem 4.1. Under hypotheses [IK3 the following formula holds 

v + = (I-2P)v~. (4.2) 

Theorem 4.2. Formula (4-*ty gives physically correct model of collision. 
Namely, denote by x(t,x,v), t G [ti , ^2] a solution such that 

x(ti,x,v)=x, x(t\, x, v) = v 

and for some x',v' and t' G (^1,^2) we have 

lim x(t,x,v') $l T x ( T i $r $i\N, x{t , x , v) G N. 
t-*r'— ' ' 

Then the solution x(t,x,v) is a continuous function of t G [^1,^2] and 
(x,v) close to (x',v'), and there exists a continuous function 

t = t(x,v), t(x,v) = t' 

such that the collision occurs at the moment r i.e. x(t(x,v),x,v) G N. 

Theorem 14.21 follows directly from the Implicit Function Theorem and 
from the fact that the function x h-> B{x) is continuous. 
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4.0.1. Proof of Theorem \4-l\ Actually we deal with vectors v € ker^4(x(7~)) 
only. By Lemma [T] one has P(ker A(x(t)) C ker A(x{t)). 

By Hypothesis [2] it follows that \v + \ = \v~\ and R(x(t)) is an isometric 
operator. 

Since the difference v + — v~ = (v^ — v7) + (v± — vj_) is perpendicular to 
ker B(x(t)) we have 




so that 

R(x(r)) |kerB(a;(T))= L 

Introduce the following space 

F(x(t)) = W{x(T))nkerA(x(T)) 

then ker A(x(t)) = F(x(t)) 0ker J B(x(r)) and R(x(t))F(x(t)) = F(x(t)). 
We finally have 

v + = Qvl + v -, Q = R( x (t))\ f{x{t)) :F(x(t))^F(x(t)). (4.3) 

Hypothesis [3] implies that Q 2 = I and consequently each eigenvalue of Q 
is either equal to 1 or to —1. 

Show that the equality Qv^_ = is possible if only = 0. Indeed, this 
equality implies v + = Qv~^_ + v jT = v~ . According to formula (|2.6p one has 
v~ € ker B(x(r)) so that vj_ = 0. Consequently Q = —I and v + = — vj_ + v7 . 
In terms of the matrix P the same is written as (|4.2p . 

The Theorem is proved. 

5. Applications 
Introduce the following notations J' = J + r 2 m, J = J + r 2 m/2. 

5.1. The Ball Rolls on the Floor and Meets the Wall. In this section 
we solve the problem we started with. 

Let tJg, or* 1 stand for velocity of the point S and for the ball's angu- 
lar velocity respectively. Superscripts + and — mark the states after the 
collision and before the collision respectively. 

In the coordinates Oxyz one has : 

^5 = v i^x + vfe y , uj^ = ujfe x + u)fe y + cofe z . 

From formula (jl.ip one has 

v f = ujfr, vf = -wfr. (5.1) 

Therefore the velocity of any point of the ball is completely defined by 
quantities v f , , uf . 
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Theorem 5.1. The following formalas hold 
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r 2 m 



■ — —v 
2J 

J' 
r,J 



rJ 



3 ) 



r 2 m 



— OJn 



2, J 



Proof. Introduce the Euler angles so that at the moment of collision one 



has cp = ifi = 0, 9 = 7r/2. Then it follows that uj 
at the moment of collision we have 



le x + ipe z — cpe v . Thus 



V-, , v , wr , uj. 



The formula T = \mv 2 s + implies G = diag(m, m, J, J, J). 



Combining formulas (15. ip and (11.21) we obtain 



.4 





1 r 



(I 









1 
1 

r 



1 



o \ 

— r 





The matrix P is calculated with the help of Lemma [TJ 
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Now the Theorem follows by direct calculation from formula (|4.2p . 
The Theorem is proved. 

5.2. The Ball is Thrown to the Floor. In this section we consider an- 
other problem with the ball. Now we have only the floor Oxy and there is 
no wall. 

Being undergone with some potential forces the ball can move in the 
half-space {z$ > r} and sometimes it can collide with the floor. 

After the ball meets the floor (zg = r) it then jumps aside. The point of 
contact C € B has the zero velocity (jl.ip . 

Introduce the configuration manifold of our system as M = M 3 x £0(3), 
where {xsiVSi^s) £ ^ 3 an d an element of 50(3) states the orientation of 
the ball. As a local coordinates in 50(3) we use the Euler angles. The floor 
is a fifth dimensional submanifold N C M which is given by the equation 
z s = r. 



Let v 



ZJ + stand for velocity of the point 5 and for ball's angular 
velocity after the collision respectively. Let C, Zd~ stand for velocity of 
the point 5 and for ball's angular velocity before the collision respectively. 
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In the coordinates Oxyz one has : 
Theorem 5.2. The following formulas hold true 
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Note that by these formulas the angular momentum about the point of 
contact C is conserved during the collision: 

m{CS,Vg] + JuJ + = m{CS,Vg] + JuJ~ . 

Due to formula (j2.9j) this is not surprise. 

Proof of the Theorem. Introduce the Euler angles so that at the moment 
of collision one has (p = tp = 0, 9 = n/2. Then it follows that uJ = 
8e x + ipe z — <pe y . Thus at the moment of collision we have 

v ± = (vf ,vf ,vf ,ujf ,uf ,ujf) T '. 

The formula T = ^mv'^ + \Joo 2 implies G = diag(m, m, m, J, J, J). From 
formula (jl.ip one obtains 



B 



10 
1 r 
,0010 




The matrix of the operator P is computed with the help of Lemma [T] 
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Now Theorem 15.21 follows from formula 
The Theorem is proved. 
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5.2.1. Nonholonomic Pendulum. Suppose that the ball moves in the stan- 
dard gravity field g = —ge z . 

Throw the ball to the floor so that 

_ _ rmv_ 

v s = —ve x — ue z: uj = — e y , u,v > 0. 

From Theorem I5J2J it follows that 



Vg = —v s , UJ = —uj . 

Thus after the ball jumped up from the floor its centre S moves along the 
same parabola just in the opposite direction. Since the angle velocity also 
changes its direction the periodic motion begins. The ball knocks the floor 
jumps up and go down along the parabola, knocks the floor at another point 
then flies along the same parabola to the initial point and so on. 

6. A Remark on Inelastic Collision 

The developed above theory allows to construct different models of in- 
elastic collision. 

For example, in the spirit of Newton's law of restitution, propose the 
following hypothesis: 

v + = — /xv J + vjT, (6-1) 

here \i £ [0, 1] is the restitution coefficient. This hypothesis is consistent 
with P~T]) . 

For plastic collision = and for super elastic one fj, = 1. In terms of 
operator P formula (16. ip has the form 

v + = (I - (1 +fi)P)v~. 

Under this hypothesis the corresponding formulas for the ball colliding 
the floor take the form 

mr 2 — uj _ Jr(l + a) _ 
v i = J, v l + J, w 2 . 

mr 2 — uj _ Jr(l + a) _ 

V 2 = J, V 2 J, W l . 

v 3 = > 

, rm(l + u) _ J — fimr 2 _ 
= J, v 2 + J, w i » 

i rm(l + a) _ J — umr 2 _ 
w 2 = J, v i + J, w 2 . 
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